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Estimate of the three-loop contribution to the QCD static potential from renormalon 

cancellation 
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It had been known that the Borel transforms of the twice quark pole mass 2mq and of the QCD 
q-q static potential V{r) have leading infrared renormalons at 6 = 1/2 such that they cancel in the 
sum. The renormalon residue of 2mq had been determined with reasonably high precision from the 
known pcrturbative coefficients of the ratio mq/niq, where rriq is the MS renormalon- free mass of 
the quark q. These values of the residues are used here to estimate the hitherto unknown part of 
the three- loop coefficient of the static potential V{r). Further, the method takes into account the 
fact that the ultrasoft energy regime contributions must be excluded from the analysis. In the bb 
quarkonium, aforementioned estimated part of the three-loop term turns out to give a contribution 
to the binding energy comparable to the nonpertubative contribution. 



This is the version v3 as it will appear in J. Phys. G. The changes in comparison to the previous 
version: the first paragraph is new; the paragraph containing Eqs. (31)-(32) is new, as is Table I; 
the last paragraph (before Acknowledgments) is new; some new references were added. Numerical 
results are unchanged. 

PACS numbers: 12.38.Bx,12.38.Cy, 12. 38. Aw 

The physics of heavy quarkonia is an interesting area of QCD because these systems are expected to be described to 
a significant degree, but not entirely, by results of perturbative QCD (pQCD). While the perturbation theory cannot 
explain the binding phenomenon, the values of the binding energies of the ground state and of some excited states of 
bb and ti are expected to be predicted to a reasonable degree by pQCD, but those of cc less well. Since masses of some 
of the heavy quarkonia [bb and cc) states are well measured, their values can be confronted with predictions based 
on pQCD, where estimates of the nonperturbative corrections should also be included. Such analyses can give us 
information on the mass of the constituents (e.g., of the heavy quark b) and predict masses of other excited states of 
quarkonium. One of the important quantities used in such analyses is the q-q static potential V(r). The perturbative 
part of V{r) is the inverted Fourier transform of the scattering amplitude for the slowly moving qq pair in the static 
limit TTtg — !■ oo. In general, however, V{r) is not reducible to a perturbative scattering amplitude that conserves all 
components of the total momentum of the two incoming and outgoing fermions on their mass-shells. In this work 
we will estimate the unknown part of the three-loop coefficient of the perturbative static potential, using the known 
property QUil^ that the leading infrared renormalon singularities cancel in the sum V{r 
pole mass of the quark. 

The perturbative expansion of the QCD q-q static potential is presently known to two loops 0,IH. In the three- 
momentum k-space, where fc = (0, k) and k = p' — p is the three- momentum transfer between the quark and antiquark 
the potential is written as [3, la, 13 



2mq, where ruq is the 
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Here, a{^) = as(/i, MS)/7r; L ~ ln(/x-^/|kp) where /i is the renormalization scale; (3j [j = 0,1,2,...) are the beta- 
coefficients appearing in the renormalization group equation (RGE) 
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We have f3o = (11 — 2n//3)/4, (3i = (102 — 38ri//3)/16. The higher f3j {j > 2) are renormalization scheme dependent, 
and are here taken in MS scheme. Here, nf = ni is the number of the hght quark flavors active in the soft/potential 
scale regime |/c^|"'^/^ ~ rUqas {uf is four for 6-6, five for t-t). The terms involving powers of L in Eq. are fixed by 
the renormalization scale independence of V. The one- and two-loop coefficients oi and a2 were obtained in j3| and 
@, respectively 
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as = 456.749 - 66.3542n/ + 1.23457n^ 
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Eus) from the soft/potential 
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Es < fJ^f < E\js ■ The existence of the IR divergent terms at ~ a in the static 



At the three- loop level [terms ^a'^ in V^(|k|) of Eq. the coefficient includes the unknown term 03 and a term with 
the infrared (IR) cutoff ^/ which cuts out the ultrasoft (US) region (|A:^|^/^ 
(S) region (|fc^|^/^ ~ m^a^ ^ Es 

Wilson loop was pointed out in Ref. 9]. The IR cutoff term 63 ln(/Ltj/|kp) in Eq. Q has 63 = 97r^/8, according 
to Ref. [inl |. and is obtained when a Green function involving color-singlet wave function of q-q is matched in two 
adjacent effective theories: nonrelativistic QCD (NRQCD) llj where the hard scales {^niq) are integrated out, and 
the potential NRQCD (pNRQCD) [l3| where the soft scales {^asTUq) are integrated out.^ 

The three-dimensional Fourier transformation of expansion gives the perturbative expansion of the static po- 
tential in the position-space 
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-a(/i)= 



^a3 + 263ln(Ai/re^^)-f /3oa2 -f -/3iai + /Sa j 2t 
I3lai + ^/3o/3i) (4^' + 7rV3) + [3l {Sf + 27r^e + 16^(3)) 



(5) 



where we use the notation £ = ln[/ircxp(7E)], with 7e the Eulcr constant (7E = 0.5772 . . .). 

The term 03 is the only part of the three-loop coefficient that has not been calculated in the literature. It is the 
value of 03 that will be estimated by using the aforementioned known property that the leading infrared renormalon 
(6 = 1/2) singularities in the sum V{r) + 2mq cancel, iriq being the pole mass of the quark. 

The 6 = 1/2 renormalon of niq tl.4j has a residue that has been determined with a reasonably high precision 
[TM IT6l|. due to a good convergence behavior of the series for the residue. First the calculation of this residue will 
be re-done, focusing on the resulting uncertainties. The perturbation expansion for the ratio S — mq/rriq — 1, where 
rfiq = rnq(^ = rnq) is the MS renormalon free mass, is known to order r^a 

^ 4 r 

= -a(^)|l + a{y) [ki -f /3oim] + a{pLf [k2 



and Lr. 
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(4/3)ki 
(4/3)k2 

ln(MV"^q)- 
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= 6.248/3o - 3.739 , 

= 23.497/3o -I- 6.248/3i + 1.019/3o - 29.94 , 
The Borel transform is thus known to order 



(2ki/3o + /3i 



(6) 
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where ri = (ki -I- /3oL,n) and r2 are the NLO and NNLO coefficients in the expansion ©, respectively. On the other 
hand, this transform can be rewritten as 



niq (1 - 26)1- 



[1 + 0(1 - 26)] + (analytic term) 



(10) 



^ In Ref. the authors calculated the N'^LO Hamilton of pNRQCD (using an estimate for of Ref. il3il ) using a method called threshold 
expansion where all the loop integrations were performed in (3 — 2e) dimensions. This generated some additional non-physical terms in 
the ultrasoft and in the soft /potential regime, but those terms canceled out in the sum from both regimes. 
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where v = /3i/(2/3q) is the power of the leading (6 = 1/2) renormalon singularity, the term in the brackets is 
independent of /i, and the last term is a function analytic in the disk |6| < 1 (the other renormalons are at 6 = 
— 1, —2, . . . ; +3/2, +2, . . ., by Refs. The /^-independent residue parameter N^, can be obtained by constructing 

the expansion of the following function in the powers of b: 



Rs{h■^Ji)^'^-{\-2bf+•'Bs{b■^JL) 

/I TT 



(11) 



This function has no pole at 6 = 1/2, only a cut at & > 1/2 which is a softer singularity. The knowledge of expansion 
Q up to results in the knowledge of expansion of Rs up to In analogy with an observation in Refs. |20|| . 

the residue parameter N,-a is 



N^^Rs{b=l/2) 



(12) 



and can thus be obtained approximately from the mentioned quadratic polynomial (truncated perturbation series - 
TPS) of Rs in b. The convergence of the TPS is good for /i ~ m,. There is some (unphysical) variation of the 
TPS predictions when /i is varied. Further, if the Pade resumniation [1/1] (6) is made, the variation with /i is 
different. The results (TPS and Pade) for N,^, as function of ii/rfiq, are given in Figs. Qla), (b) for Uf — 4,5, 
respectively. The points of zero /i-sensitivity (principle of minimal sensitivity - PMS) give us in the case = 4 
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FIG. 1: The residue parameter Nm determined from (N'^LO) TPS and Pade [1/1] of Rs(b) (6=1/2), as function of the renormahzation 
scale parameter x = fi/rrlq, for rif = 4 (a) and nj = 5 (b). 

predictions Nm — 0.563,0.546, for TPS and Pade, respectively, and in the case nf = 5 predictions Nm = 0.539,0.526. 
As central values we take the arithmetic average between the TPS and Pade values at respective zero /i-sensitivities: 
Nm = 0.555(n/ = 4);0.533(n/ = 5). We could also include in the analysis the Pade predictions at such ^ where 
bp^i^ = —1, i.e., the theoretical nearest-to-origin pole of i?s(6). However, interestingly, such /i's are close to the 
Pade-PMS /i's {p/rrTq 0.85 — 0.87 when Uf = 4; n/ruq « 0.90 — 0.95 when Uf = 5), and the predictions for Nm 
are virtually the same. The TPS predictions for Nm start falling down fast when ix/rriq decreases below the value 
one. However, [ijrriq < 1 when the Pade approximant gives the value 6poic = —1 ifJ-fT^q = 0.85,0.90 for Uf = 4,5, 
respectively). The absolute value of the uncertainty will be estimated as the deviation of the TPS prediction of Nm 
at such low /i from the aforementioned central values (0.555,0.533). The resulting estimates are thus 

Nm{nf = 4:) = 0.555 ±0.020, (13) 
Nm{nf = 5) = 0.533±0.020. (14) 

On the other hand, the TPS predictions for Nm at ^/rfiq = 1 agree exactly with the values obtained by Pineda [Tsf 
{Nm = 0.5523,0.5235, for Uf = 4,5, respectively), and are well within the limits H13() - (|14|l . Further, in Ref. "T^ the 
values Nm ~ 0.557 ± 0.008, 0.530 ± 0.008 for = 4, 5 were obtained using a method involving a conformal mapping, 
and they are close to the values H13|l -H14 () here, but the estimated uncertainties here are larger. 

A somewhat analogous procedure is now applied to the three-loop expression |(SJ) for the static potential: Requiring 
the reproduction of the values of Nm lfO)l - l(T^ by the static potential lO will allow us to estimate the values of the 
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unknown three-loop coefficient expression 03 appearing in expansion (jSJ 

53^| + ^ln(A./re^-) (15) 

as a function of the renormalization scale parameter X = fxr (it is to be recalled that 63 = 97r^/8). Defining the 
dimensionless quantity F{t) = (— 3/(47r))rF(r), the TPS (0) gives us the Borel transform of F{r) up to order 

where Vj are the coefficients at powers of a{^) in expansion jsj (vi = ai/4 + 2/3o^, etc.). This function has the 
renormalons at 6 = 1/2,3/2,5/2, etc. [23. It can be written as 

3 

Bpib; fir) ^-fir _ [1 + 0{1 - 2b)] + (analytic term) , (17) 

where 1/ = /3i/(2/3q), the expression in the brackets is /i- independent, and the last term is analytic for |&| < 3/2. 
The parameter Nm here is the same as in expression (|l()(l to ensure the leading renormalon cancellation in the Borel 
transform of the sum (y{r) + 2mq). The function 

Rp{b; ^r) = ~{l- 2b)^+''BF{b- fir) , (18) 
3 iir 

is then less singular at b = 1/2 (cut instead of pole), and its evaluation at 6 = 1/2 should give us 

RF{b=ll2-nr)=Nm . (19) 

The coefficients vi and V2 appearing in expansions jSJ and (|16|l are explicitly known, but W3 contains the undetermined 
three-loop parameter 03 (|15|l a^]. This implies that the coefficients in the expansion of Rpib) are known explicitly 
up to terms but the coefficient at 6^ contains 03. Thus, any resummation of the N'^LO TPS of Rpib) (i.e., the TPS 
including the &'^-term) at 6 = 1/2, and requiring the identity p9|l . with Nm values H13|l -H14 |l . gives us in principle an 
(approximate) value of the unknown three- loop coefficient 03. 

From the purely perturbative point of view, the procedure indicates inconsistency because we are matching the 
residue parameter obtained on the basis of the N^LO TPS for Rs{b) (|ll|l with the residue parameter obtained on 
the basis of the N^LO TPS for Rpib) (|18|l . However, in practical terms, there appears to be no inconsistency. The 
reason lies in the considerably better convergence behavior of the N^LO TPS of Rs{b) than the N^LO TPS of Rpib) 
at 6 = 1/2. The smallness of the term in Rs{b) does not seem to be accidental, because this term is small for any 
relevant value of [uf < 5). Furthermore, for these reasons, the procedure of estimating 03 is indeed rather stable, 
because the N'^LO term (^6"^) in Rp^b) now plays an important role in ensuring numerically the 6=1/2 renormalon 
cancellation condition H19I) . 

On the other hand, since RF(b) is a significantly less singular function than Bp{b) and at the same time the TPS 
of Rp{b) shows only slow convergence (at b = 1/2), it is very reasonable to apply some form of resummation beyond 
the simple TPS evaluation - e.g., Pade resummation. The most stable results for 03, when fi is varied, are obtained 
by using the [2/1] (6) Pade of (the N^LO TPS of) Rp{b). The resulting values of the parameter 03, as function of 
the renormalization scale parameter X = fj,r and requiring the central values Nm of Eqs. (|13|I " H14|I . are presented in 
Figs. |2Ia),(b), for the cases Uf = 4,5, respectively. When Pade [2/1] resummation is apphed to Rp{b), the PMS 
points do exist, with the values 0.3 = 77.9, 45.2 at ApMS ~ 0.56, 0.60, when Uf = 4, 5, respectively. If the values of Nm 
increase, so does 03. The final result is 

53(ri/=4) = 77.9l^l'l{SNm)±U.2{Sti) , (20) 
53(n/ = 5) = A5.2 t llidNm) ± 9.1(6^1). (21) 

where the first uncertainty is due to the uncertainty ±0.020 in the values of Nm (|13l) ~ (|14|l (Pade [2/1] was evaluated at 
the unchanged aforementioned values of A), and the second uncertainty is due to the renormalization scale dependence. 
The latter uncertainty was evaluated as plus / minus the deviation of 03 from the mean value when A = /ir is increased 
by factor two from its PMS value. 

The results l)2(J |l - (|21() for the three-loop parameter 03 (|15|) can be extended further, in order to estimate the three- 
loop coefficient 03. This will involve further uncertainties of the S-US factorization scale /i/ and of the typical distance 
r between the quarks. As emphasized in Ref. 10], the considered q-q singlet static potential is not the quantity defined 
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FIG. 2: The three-loop parameter as determined from N^LO TPS and Fade [2/1] of Rpib) (6 = 1/2) and requiring the central values 
1131 — il4i for the residue parameter Nm, as function of the renormalization scale parameter A = ^r, for nf = 4 (a) and nj = 5 (b). 



via the vacuum expectation value of a static Wilson loop which includes US contributions, but rather the relevant 
object for the dynamics of the q-q pairs with large but finite mass. The kinetic and the binding energies of the q-q 
system are US energies. Therefore, the interactions of the US gluons with the q-q system are sensitive to the energies 
of that system and should be excluded from the static potential 10]. The exclusion of these interactions leads to the 
term oc \n{fifre'^^) at in the potential Hence the value of the parameter 03 will depend on the values of the 
scales and r. The dependence of V{r) on the IR cutoff ^/ of the soft scale regime, and the consequent dependence 
of the values of Nm (|19|l and as [cf. Eq. H15|) ] on /i/, does not imply that US contributions affect them, but rather 
that the US contributions to those quantities are cut out, zero. 

Since the scales /^/ and r appear in logarithms, a natural choice for /i/ would be /i/ = {EgEus 
S scale is of the order of the three-momentum transfer |k| in the q-q ground state |1): Eg ^ mqas 
scale is i?us ~ l^qql ~ mqal- Therefore, the factorization scale can be estimated as 



1/2. The typical 
The typical US 



(EsE^s) 



1/2 _ 



(22) 



where k ^ 1 and we take /is ~ M (~ -^s)- The typical distance r in the combination 03 (|15|l can be estimated in 
the following way. First the expectation value of the static potential term proportional to 03 is calculated in the q-q 
ground state |1} 



{l\-a%f,)Mr)\l) 
r 



1 



1 



03 Qtt^ 



In (/i/aB(Ai)e/2) 



03 Stt- /3e k 1/2/ N 



(23) 
(24) 



Here, asifJ-) — 3/(2mgas(/i)) is the modified Bohr radius |22j | appearing in the qq ground state wavefunction, and we 
denoted 



K= — -— = I + a[n] 



a(/i)2 



^02 + (ai/3o + 2/3i) 7E + Po (47! + ^Vs) 



(25) 



In Eq. (|24|l . the expression H22() for /i/ was inserted. Eqs. (|23|I - H24|) suggest that we can estimate the value of 03 from 
the value of 03 via the relation 



Ia3«a3-^ln(2.04 ^al^'if,. 



(26) 



Now we estimate the values of k and k for the b-b and t-i system. 

Since /i^ ^ mqas{fis), we have for the b-b system (m^ « 5. GeV, Uf = ni = 4), the values fig ~ 1-2 GeV and 
as{fJ-s) ~ 0.3 (as suggested by resummations of the semihadronic T-decay width [23). The series 1^3)1 is strongly 
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TABLE L The separate contributions of the perturbative potential V{r) @ for two typical radia r of the bb system (uf —4). 
The as-contribution is given for the central value (1291 of 03. The values are given for the renormalization scale p — 2/r, and 
for fi = 1/r in parentheses. Included are also values of the logarithmic confining potential of Ref. 2(3]|. All energies are in MeV; 
radia are in GeV^^. Other details are given in the text. 



r 


Konf 


V{as) 


V(LO) 


V(NLO) 


^(N^LO) 


V(N^LO) 




0.50 
0.75 


+19 
+43 


-23( -74) 
-29(-126) 


-629(-842) 
-491(-710) 


-319(-328) 
-292(-350) 


-281(-374) 
-301(-505) 


-326(-553) 
-409(-944) 


0.236(0.316) 
0.276(0.399) 



divergent then (1 + 0.37 + 0.40 +•••), and we can estimate /5(n/ = 4) w 1.25 ± 0.10 by ap plyi ng the resummation 
method of ,2j|.^ The resummed binding energy in the ground state is E'jj « —0.4 GeV l25l|. i.e., i?us ~ l-E-febl ~ 
0.9 iTibasins)'^ ■ On the other hand, Es is approximately equal to the expectation value of the three-momentum in the 
ground state 

i^s-(l||p||l) = ^^ = g^™,a.(.). (27) 

Therefore, Eq. gives k « 0.7^^/^ (« 0.8) when n/=4. 

In the case of the t-t system (m^ « 175 GeV, n-f = rii ~ 5), we have fis ~ 30 GeV, and as{ns) ~ 0.14. The series 
is less divergent (1 + 0.15 + 0.07 + •■• ), and we estimate with the aforementioned resummation K{nf — 5) « 1.30 + 0.10. 
The resummed binding energy in the ground state is Eu « —3. GeV 0,113, i.e., E\js « \Ett\ ~ 0.9 mtas{fJ.s)^- 
Eqs. ^ and then give k « O.Ik^^'^ (« 0.8) when Uf^d. 

Therefore, the approximate relation (|26|l can then be rewritten explicitly for the b-b and t-t system as 

^«3«a3-^ln(l.4^ay2(^,)j , (28) 

where, as mentioned, we take the central values: as{fJ-s) ~ 0.3 and k « 1.25 for rt/=4; as{fJ-s) ~ 0.14 and k w 1.3 for 
nf = b. This then gives 

^ a3(n/-4) - 86.3 t,l\m,,,) ± li.2{5^i) t I'.HSfif) , (29) 

^ a3(n/ = 5) = 62.5 t ^JSN.^) ± 9.1(6 ^i) ± 15.4(<5m/) , (30) 

where the last uncertainty is from the uncertainty in the argument of the logarithm in Eq. (|28|) involving the estimates 
of the soft and ultrasoft scales. These uncertainties are large, because the central values in the logarithm of Eq. H28|) 
are only suggestive. The value ±15.4 for these uncertainties was obtained by allowing the overall factor 2 or 1/2 in 
the logarithm in Eq. H28|l . On the other hand, the smaller downward uncertainty —8.4 in the case n/ =4 is due to the 
problematically small hierarchy between the soft and ultrasoft scales in this case, and the related minimal requirement 
that the factorization scale fif should be below the typical soft scale ~ |k|. Stated otherwise, the logarithmic term 
in Eq. (|15|l should be negative, or equivalently, the IR cutoff term proportional to 63 in expansion |SJ| should give a 
positive contribution to V{r) (i.e., fif < 0.56/r). 

We now investigate how numerically important the obtained a3-terms are with respect to other terms in the binding 
energy of the quarkonium. For example, in the bb ground state T(IS'), the typical value of the interquark distance is 
r l/{mbas{fJ,s)), i-e., r « 0.50-0.75 GeV~^. TablelJshows the values of the pure as-term of V{r) for such r, for the 
central 03 value (|29l) . as well as the values of the entire leading order (LO; ~ o/r), NLO (~ /r), N^LO (~ a'^ /r), 
and N^LO a'^/r) parts of V{r) of Eq. ||SJ). These values are given for two choices of the renormalization scale: 
fi = 2/r, 1/r. The values of the a^-term as given in the Table should be regarded with some caution though, because 
they strongly depend on the choice of fi and because the perturbative series for V{r) is asymptotically divergent. This 
can be seen from the values of different orders of V{r) given in the Table. Borel-related resummation methods for the 
binding energy £'(,5 |25| give for the a3-term contribution (—16 + 4) MeV, i.e., £^(,5(03) — £'(,5(0) « (—16 + 4) MeV, when 



^ Note that the series 1251 represents the quantity (—3/{4iT))rV{r) of Eq. J5J at fir = 1. The principal value (PV) prescription was taken 
for the Borel integration over the b = 1/2 renormalon singularity. 
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the values of 03 are taken from Eq. 129|l ^ and jj, ^ Eg. In addition, we display in TableUthe values of the logarithmic 
confining potential V^onf ('') — (CpaC/'K)Vo(Cr) of Ref. (26i |. based on the renormalization-group- improved light-front 
Hamiltonian formalism of Refs. I27ll28l 



4 

K;onf(?') ^ —aC 

Stt 



21ni?- 2Ci(i?) + ^Si{R) - -^(1-cosi?) + -| sin i? - 5 + 27e 



(31) 



where: R = Cr; >C is a cutoff scale parameter of the framework; a is an effective as- Ci and Si are the usual cosine and 
sine integrals: Ci(z) = — f [cos{t)/t]dt; Si(z) = [sm{t)/t]dt. The fitted values of the parameters in that potential 
are also taken from Ref. |2l|: c=l [c = (4/3)(m6a//:)], mf, = 4.8 GeV, a = 0.4 £ = 2.56 GeV). The nonperturbative 
contribution to the binding energy Ef^i can alternatively be estimated as coming from the giuonic condensate |29j| 

£;,5(7/s)("P) « 7716^2 gl (^^nTbasifius)^ (a(/^«.)Gp.G^") « (50 ± 35) MeV , (32) 

where Tni, — 4.2 GeV was taken, and the ultrasoft energy ii^s ~ mf,a^ < 1 GeV was taken equal to /i « 1.5-2.0 GeV in 
order to be able to determine as{fJ,us) still perturbatively [as(A*) ~ 0.30-0.35]. Thus, the confining energies of Table 
U and of Eq. H32() consistently give nonperturbative contributions to Ef^i roughly in the range 20-50 MeV, compared 
to (—16 ± 4) MeV of the oa-term. We see that the aa-term in bb gives a nonnegligible contribution in comparison 
to the nonperturbative contributions to the binding energy (and thus to the meson mass). On the other hand, for 
the lighter charmonium the nonperturbative effects are expected to be more important, while for the much heavier 
toponium less important than the aa-term contributions. 

The contributions to Ef,^ from the "hard" (~ m&) and "soft" (~ rubas) energy regimes give about —300 MeV, and 
those from the ultrasoft regime (~ m^a^) about —150 MeV |2a|. By the virial theorem, this implies for the kinetic 
energy an estimate of 150-250 MeV. 

We now compare the results (|20ll - (|21|) and (|29|I " H30|I with those in the literature. The authors of Ref. J^] have 
estimated the three-loop coefficient of the static potential V^(|k|) with a Pade-method based on the known one- and 
two-loop coefficients. However, they did not have the IR cutoff term in the three-loop coefficient. Therefore, there is 
no direct correspondence with the results presented here, but in V{r) their coefficient cq formally may correspond to 
03 of Eq. fTB|l . or equivalently, to 03/4^ if fif = exp(— 7E)/r. Their estimated values were cq = 97.5, 60.1 for n/ = 4, 5, 
which are higher than the values (|20ll -121 |l . The estimates of T. Lee [l^, also without the effects of the IR cutoff, and 
based on a Borel transform method j3Qj, give 03/4'^ = 59±81(n/ = 4); 34±63(n/ = 5). Further, Pineda 0| estimated 
the three-loop coefficient of V{r) from renormalon-dominated large order behavior of the coefficients for ^ = 1/r but 
without consideration of the IR cutoff effects; if his results are to be interpreted with the value of the factorization 
scale fif = 1/r {fif = v^s in the notation of 15]), as suggested by him, then this would give, in the notation l|15|) 
of the present paper, the values 03/4"^ = 59.6,24.3 for Uf — 4,5. On the other hand, in terms of 03 of Eq. I|15|) . his 
results would imply 03 — 72.4,37.1 for Uf = 4,5, which are not very far from the values (|20|l -121 |) . Pineda, in the 
framework of his method, did account for the leading IR renormalon cancellation in {V{r) + 2mq). 

The main results and conclusions of the present work are the following: The unknown part 03 of the three-loop 
coefficient of the static potential was estimated by using the known property that the leading infrared renormalon 
singularities cancel in the sum V{r) + 2mq. Further, the presented method takes into account the fact that the 
contributions of the ultrasoft energy regime should be excluded from the analysis. The obtained estimated values of 
03 for the bb and tt quarkonia are given in Eqs. and (EOJ, respectively. In the bb system, the obtained value of 
the three-loop parameter 03 leads to a decrease of the binding energy by about 10-20 MeV. This is smaller than, but 
still comparable to, the increase of the binding energy by about 20-50 MeV by nonperturbative effects. The similar 
value (|30|l of 03 for the much heavier tt system is expected to influence the binding energies to a similar degree as in 
bb, but the nonperturbative effects are expected to be much weaker. For the lighter cc system the nonperturbative 
effects are expected to dominate over the three-loop effects. 



The uncertainties were added in quadrature: 03/4'^ Ri 86 ± 23. 
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